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Abstract. The decay J/tp — > copp is studied using a 5.8 x 10 7 J/tp event sample accumulated with the 
BES II detector at the Beijing electron-positron collider. The decay branching fraction is measured to 
be B(J/ip — > copp) = (9.8 ± 0.3 ± 1.4) x 10~ 4 . No significant enhancement near the pp mass threshold 
is observed, and an upper limit of B{J/ip -» uA(1860))B(X(1860) -> pp)< 1.5 x 10~ 5 is determined at 
the 95% confidence level, where X(1860) designates the near-threshold enhancement seen in the pp mass 
spectrum in J/ip — > ypp decays. 
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1 Introduction 

Decays of the J/tp meson are regarded as being well suited 
for searches for new types of hadrons and for systematic 
studies of light hadron spectroscopy. Recently, a number of 
new structures have been observed in J/tp decays. These 
include strong near-threshold mass enhancements in the 
pp invariant mass spectrum from J/ip —> jpp decays [I], 
the pA and K~ A threshold enhancements in the pA and 
K~A mass spectra in J/tp — > pK~A decays [2], the ui(f> 
resonance in the Lu<p mass spectrum in the double-OZI 
suppressed decay J/tp — > juj(f> [3], and a new resonance, 
the (1835), in J/tp -> 7^+^-77' decays g]. 

The enhancement X(1860) in J/tp — > -fpp can be fitted 
with an S- or P-wave Breit-Wigner (BW) resonance func- 
tion. In the case of the S-w&ve fit, the mass is 18591^0^25 
MeV/c 2 and the width is smaller than 30 MeV/c 2 at the 
90% confidence level (C.L.). It is of interest to note that 
a corresponding mass threshold enhancement is not ob- 
served in either pp cross section measurements or in B- 
meson decays [5]- 

This surprising experimental observation has stimu- 
lated a number of theoretical interpretations. Some have 
suggested that it is a, pp bound state (baryonium) [6,7,8, 
[9lfTfj], Others suggest that the enhancement is primarily 
due to final state interactions (FSI) between the proton 
and antiproton [111112]. 

The CLEO Collaboration published results on the ra- 
diative decay of the T(1S) to the pp system [13], where 
no pp threshold enhancement is observed and the up- 
per limit of the branching fraction is set at B(T(1S) — ► 
7 X(1860)) J B(X(1860) -> pp) < 5 x 10" 7 at 90% C.L.. 
This enhancement is not observed in BES2 tp(2S) — > jpp 
data either [T3] and the upper limit is set at B{tp(2S) — > 
7X(1860))B(X(1860) pp) < 5.4 x 10~ 6 at 90% C.L.. 

The investigation of the near-threshold pp invariant 
mass spectrum in other J/tp decay modes will be helpful 
in understanding the nature of the observed new struc- 
tures and in clarifying the role of pp FSI effects. If the 
enhancement seen in J/tp — > jpp is from FSI, it should 
also be observed in other decays, such as J/tp — ► upp, 
which motivated our study of this channel. In this paper, 
we present results from an analysis of J/tp — > 7r + 7r~7r°pp 
using a sample of 5.8 x 10 7 J/tp decays recorded by the 
BESII detector at the Beijing Electron-Positron Collider 
(BEPC). 

BES is a conventional solenoidal magnetic detector 
that is described in detail in Ref. [15]. BESII is the up- 
graded version of the BES detector [IB]. A twelve-layer 
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Vertex Chamber (VC) surrounds a beryllium beam pipe 
and provides track and trigger information. A forty-layer 
main drift chamber (MDC) located just outside the VC 
provides measurements of charged particle trajectories over 
85% of the total solid angle; it also provides ionization en- 
ergy loss {dE/dx) measurements that are used for particle 
identification (PID). A momentum resolution of a p /p = 
1.78%a/1 + p 2 (p in GcV/c) and a dE/dx resolution of 
^8% are obtained. An array of 48 scintillation counters 
surrounding the MDC measures the time of flight (TOF) 
of charged particles with a resolution of about 200 ps for 
hadrons. Outside of the TOF counters is a 12 radiation 
length, lead-gas barrel shower counter (BSC), that oper- 
ates in self quenching streamer mode and measures the en- 
ergies and positions of electrons and photons over 80% of 
the total solid angle with resolutions of <Je/E = 0.21/^/E 
(E in GeV/c 2 ), er^ = 7.9 mrad, and a z = 2.3 cm. External 
to a solenoidal coil, which provides a 0.4 T magnetic field 
over the tracking volume, is an iron flux return that is 
instrumented with three double-layer muon counters that 
identify muons with momentum greater than 500 MeV/c. 

Monte-Carlo simulation is used to determine the mass 
resolution and detection efficiency, as well as to estimate 
the contributions from background processes. In this anal- 
ysis, a GEANT3-based Monte-Carlo program (SIMBES), 
with a detailed simulation of the detector performance, is 
used. As described in detail in Ref. [17], the consistency 
between data and Monte-Carlo has been validated using 
many physics channels from both J/tp and tp(2S) decays. 



2 Analysis of J/ip — > upp, to — > n + it n° 

For candidate J/tp — > 7r + ir~ ir c 'pp events, we require four 
well reconstructed charged tracks with net charge zero in 
the MDC and at least two isolated photons in the BSC. 
Each charged track is required to be well fitted to a he- 
lix, be within the polar angle region | cos f?| < 0.8, have 
a transverse momentum larger than 70 MeV/c, and have 
a point of closest approach of the track to the beam axis 
that is within 2 cm of the beam axis and within 20 cm 
from the center of the interaction region along the beam 
line. For each track, the TOF and dE/dx information is 
combined to form a particle identification confidence level 
for the 7T, K and p hypotheses; the particle type with the 
highest confidence level is assigned to each track. The four 
charged tracks are required to consist of an unambiguously 
identified p, p, Tt + and 7r~ combination. An isolated neu- 
tral cluster is considered as a photon candidate when the 
angle between the nearest charged track and the cluster 
is greater than 5°, the angle between the p track and the 
cluster is greater than 25° [TH], the first hit is in the be- 
ginning of six radiation lengths of the BSC, the difference 
between the angle of the cluster development direction in 
the BSC and the photon emission direction is less than 
30°, and the energy deposited in the shower counter is 
greater than 50 MeV. A four-constraint kinematic fit is 
performed to the hypothesis J/tp — > pj37r + 7r~77, and, in 
the cases where the number of photon candidates exceeds 
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two, the combination with the smallest % 2 
is selected. We further require that Xpp„+„- 

Figure [T] shows the 77 invariant mass of the events 
which survive the above-listed criteria, where a distinct 
7T° — * 77 signal is evident. Candidate tt° mesons are se- 
lected by requiring |M 77 — m w o| < 0.04 GeV/c 2 . After this 
selection, a total of 15260 events is retained. The 7r + 7r~7r° 
invariant mass spectrum for these events is shown as data 
points with error bars in Fig. [21 where prominent lo and r) 
signals are observed. 




0.05 0.1 0.15 0.2 0.25 0.3 

M(yy) (GeV/c 2 ) 

Fig. 1. The M 77 distribution for J/ip — » jjtt + iy~pp candidate 
events. 




M(jrW) [GeVIc? 



Fig. 2. The M % + n - n o distribution for J/tp — > tt + -k~tt°pp can- 
didate events. The dots with error bars are data. The solid 
histogram is the background estimated from Monte-Carlo sim- 
ulation, normalized according to the PDG branching fractions. 
The solid curve is the result of a fit described in the text. The 
dashed curve is the background polynomial. 



The backgrounds in the selected event sample are stud- 
ied with Monte-Carlo simulations. We generated J/tp — > 
ppir + Tr~ir° decays as well as a variety of processes that 
are potential sources of background^ J/tp — > ppr)'(r)' — ► 
TT + ir-ri);ppr)'(ri' -> p°j); ppir+ir-; AAir"; £°£ a ; E(l385)-£+; 
?7 C 7; A++A — ; jppTr + ir-; A ++ pTT~; AS^tt+ (+ c.c); E°tt°A; 
r(1385)°r°; A ++ A--Tr a ; and in proportion to 

the branching fractions listed in the Particle Data Group 
(PDG) Tables |19j . The main background sources are found 
to be the decays J/tp -> AS^tt + (+ c.c.) and A++A 7r°. 
The 7r + 7T _ 7r° invariant mass spectrum for background events 
that survive the selection criteria is shown as a solid his- 
togram in Fig. [2 here no signal for Lopp is evident. 

The branching fraction for J/tp — > ojpp is computed 
using the relation 



B{J/tp — > ujpp) — 



N, 



obs 



N 



TT+TT-TT ) • B(TT° 



■ 77)' 



Here, N obs is the number of observed events; Nj/^ is 
the number of J/tp events, (57.7 ± 2.6) x 10 6 [20 ; e is the 
Monte-Carlo determined detection efficiency; and B(oj — ► 



~7T°) and B(ir — > 77) are the u> 
77 branching fractions. 



and 



The 7T- 



invariant mass spectrum shown in Fig. [5] 



is fitted using an unbinned maximum likelihood fit with 
resolution broadened BW functions to represent the lo and 
r) signal peaks. The mass resolutions are obtained from 
Monte-Carlo simulation to be 12 MeV/c 2 for the u> and 
and 14 MeV/c 2 for the ij. The masses and widths of the 
lo and r] are fixed at their PDG values [19] . A 4th-order 
Chcbychev polynomial is used to describe the background. 
The fit gives an to signal yield of 2449±69 events. The de- 
tection efficiency from a uniform-phase-space Monte-Carlo 
simulation of J/tp — > ujpp [lo — > 7r + 7r~7r , 7r° — * 77) is 
4.9 ± 0.1)%. The branching fraction is determined to be: 

B(J/tp -> ujpp) = (9.8 ± 0.3) x 10~ 4 , 

where the error is statistical only. 

We use this sample with \M 7r + n - 7T a — 0.783| < 0.03 
GeV/c 2 to study the near-threshold region of the pp in- 
variant mass spectrum. Figure [3] shows a Dalitz plot for 
the selected J/tp — > ujpp candidates, where no obvious 
structure is observed although it is not a uniform distri- 
bution. Figure [4] shows the threshold behavior of the pp 
invariant mass distribution. The dotted curve in the figure 
indicates how the acceptance varies with invariant mass. 

The backgrounds in the pp threshold region mainly 



come from the decays of J/tp — > AS n + (+ 



c.c. 



and 



A++A- 



The M(pp) dependence of this background 



can be modeled by appropriately scaled data from the u> 



sidebands (0.663 GeV/c 2 < M d 



< 0.723 GeV/c 2 and 



0.843 GeV/c 2 < M„ 



<0.903 GeV/c 2 ). 



The contributions of sideband and non-resonant ujpp 
events can be well described by a function of the form 



f(S) = N(5i + Oi^3 + a 2 S^) 



with S — M p p 



2m v 
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Fig. 3. The Dalitz plot for J/ip — ► copp candidate events. 

In Fig. [31 no significant excess over the background 
plus non- resonant terms is evident. A Bayesian approach [TU] 
is employed to extract the upper limit on the branching 
fraction of J/ip — > ujX (1860). An acceptance-weighted S- 
wave BW function 

(2;+i) fc 3 

BW M« (M*-M3r-M3r> - e M 

is used to represent the low-mass enhancement. Here, r 
is a constant width, q is the momentum of proton in the 
pp rest frame, I is the relative orbital angular momen- 
tum of p and p, k is the momentum of u>, and e(M) is 
the detection efficiency obtained from Monte-Carlo sim- 
ulation. The mass and width of the BW signal function 
are fixed to 1860 MeV/c 2 and 30 MeV/c 2 , respectively. 
The contributions of background and non-resonant copp 
events are presented by the function form f(5), where 
the parameters a\ and are allowed to float. As shown 
in Fig. [U the solid curve is the fit of the M p p - 2m p 
with the BW signal function and f(S) function described 
above. Using the Bayesian method, the 95% C.L. upper 
limit on the number of observed signal events is 29. Since 
the J PC of X(1860) is unknown, we use simulated events 
distributed uniformly in phase space to determine a de- 
tection efficiency of J/ip -> w AT (1860) (X(1860) -> pp, 
uj -> tt+tt-tt , tt° -> 77) of (4.7 ± 0.1)%. The upper limit 
of the branching fraction, without considering the system- 
atic errors, is then: 

B(J/ip -> wX(1860)) • B(X(1860) -> pp)) 



Nj/tp ■ e ■ B(ui -> tt+tt-tt ) • B(tt° -> 77) 



3 Systematic errors 

The systematic errors on the branching fractions are mainly 
due to uncertainties in the MDC tracking, kinematic fit- 




Fig. 4. The M p p — 2m p distribution for J/ip —* copp candidate 
events. The dots with error bars are data. The solid curve is 
the result of fit described in the text. The dashed curve is the 
function used to represent the background plus non-resonant 
Lopp events. The dotted curve indicates how the acceptance 
varies with pp invariant mass. 

ting, particle identification (PID), photon detection, back- 
ground estimation, the model used to describe hadronic 
interactions in the material of the detector, and the un- 
certainty of the total number of J/ip decays in the data 
sample. 

The systematic error associated with the tracking ef- 
ficiency has been carefully studied [T7]. The difference of 
the tracking efficiencies between data and Monte-Carlo is 
2% per charged track; an 8% contribution to the system- 
atic error associated with the efficiency for detecting the 
four-track final state is assigned. In Ref. [P71I21] . the effi- 
ciencies for charged particle identification and photon de- 
tection are analyzed in detail. The systematic errors from 
PID and photon detection are 2% per proton (antipro- 
ton), 1% per pion and 2% per photon. In this analysis, 
with four charged tracks and two isolated photons; 6% is 
taken as the systematic error due to PID and 4% due to 
photon detection. The uncertainty due to kinematic fitting 
is studied using a number of exclusive J/ip and ip(2S) de- 
cay channels that are cleanly isolated without a kinematic 
fit [2"2l23] . It is found that the Monte-Carlo simulates the 
kinematic fit efficiency at the 5% or less level of uncer- 
tainty for almost all channels tested. Therefore, we take 
5% as the systematic error due to the kinematic fit. 

The background uncertainties come from the uncer- 
tainty of the background shape. For the branching frac- 
tion measurement of J/ip — * copp, changing the order of 
the polynomial background causes an uncertainty in the 
number of background events. For the upper limit deter- 
mination of J/tp — * ujX (1860), the uncertainty of back- 
ground shape can be determined by the fitting results 
with the background shape fixed to the function form 
f(S), derived from fitting the scaled ui sideband data plus 
phase-space generated ujpp MC events. Respectively, 5% 
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and 10% are taken as the systematic errors due to the 
background uncertainties in the branching fraction mea- 
surement of J/tp — > ujpp and the upper limit determina- 
tion of J/tp -> wl( 1860). 

Different simulation models for the hadronic interac- 
tions in the material of the detector (GCALOR/FLUKA) [H 
[25] give different efficiencies. Respectively, 4.8% and 11.4% 

are taken as the systemic errors due to the different hadronic 5 Acknowledgments 

models in the branching fraction measurement of J/tp — > 
wpp and the upper limit determination of J/tp — > wAT(1860) 
In addition, if the J p of X(1860) is 0~ , the angular distri- 
bution of the u> would be l + cos 2 8. A Monte-Carlo sample 
generated with the to produced with a 1 + cos 2 9 distribu- 
tion and a uniform distribution for the X(1860) decay into 
pp results in an 8.5% reduction in detection efficiency. This 
difference is taken as the systematic error associated with 
the production model. 

The branching fractions of u> — ► 7r + 7r~7r and tt° — > 77 
are taken from the PDG tables. The errors of the interme- 
diate decay branching fractions, as well as the uncertainty 
of the number of J/tp events [20] also result in the system- 
atic errors in the measurements. 

The systematic errors from the different sources are 
listed in Table[T] The total systematic errors for the branch- 
ing fractions are obtained by adding up all the systematic 
sources in quadrature. 



is determined. The absence of the enhancement AT (1860) 
in J/tp — ► wpp, T(1S) — > 7pp and tp(2S) — ► jpp also 
indicates its similar production property to that of 77' [261 
[27] . i.e., AT(1860) is only largely produced in J/tp radiative 
decays. 
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Table 1. Systematic error sources and contributions (%). 





~B(J/ip — > Lupp) 


Upper Limit 


Tracking efficiency 
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4 Summary 
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B(J/tP -> tupp) = (9.8 ± 0.3 ± 1.4) x 10~ 4 . 
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